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ABSTRACT 

In this paper we derive the mass function of seed black holes that result from the 
central mass concentrated via disc accretion in collapsed haloes at redshift z w 15. 
Using standard arguments including stability, we show that these pre-galactic discs 
can assemble a significant mass concentration in the inner regions, providing fuel 
for the formation and initial growth of super-massive black holes. Assuming that 
these mass concentrations do result in central seed black holes, we determine the 
mass distribution of these seeds as a function of key halo properties. The seed mass 
distribution determined here turns out to be asymmetric and skewed to higher masses. 
Starting with these initial seeds, building up to 10^ solar masses by z = 6 to power 
the bright quasars is not a problem in the standard LCDM cosmogony. These seed 
black holes in gas rich environments are likely to grow into the supermassive black 
holes at later times via mergers and accretion. Gas accretion onto these seeds at high 
redshift will produce miniquasars that likely play an important role in the reionization 
of the Universe. Some of these seed black holes on the other hand could be wandering 
in galaxy haloes as a consequence of frequent mergers, powering the off-nuclear ultra- 
luminous X-ray sources detected in nearby galaxies. 

Key words: accretion, accretion discs - black hole physics - galaxies: formation - 
cosmology: theory ~ instabilities ~ hydrodynamics 



1 INTRODUCTION 



It is becoming increasingly clear that supermasssive 
black holes assemble in galactic nuclei, and that their 
;rowth is in tandem with that of the stellar component 
Merritt fc Ferrarese|[200ll : iTremaine et al.ll2002l ). Optically 
bright quasars powered by accretion onto black holes are 
now detected out to redshifts of z > 6 w hen the Universe 
was b arely 7% of its current age (|Fan et al . 2004; F an et al.l 
l2006f ) . The luminosities of these high redshift quasars imply 
black hole masses Mbu > 10^ Mq. Assembling these large 
black hole masses by this early epoch starting from remnants 
of the first generation of metal free stars has been a challenge 
for models. Models that describe the growth and accretion 
history of supermassive black holes com mence with high- 
redsh i ft seeds derived from Pop-Ill stars |Haiman fc Loebl 
1 19981 : IVolonteri et"all l2005l : IWvithe fc Loebl l2005h . Some 
suggestions to accomplish this invoke sup er-Eddington ac- 
cretio n rates for brief periods of time (iVolonteri fc Ree3 
I2OO5I ). Alternatively, it has been sugge sted that the for- 
mation of more massive seeds ab-initio (|Eisenstein fc Loebl 
I1995I : iKoushiappas et al] |2004|'). possibly through t he col- 



also alleviate the problem of building up supermassive black 
hole masses to the required values by z = 6. We explore this 
possibility further in this paper. 

In a previous paper (jLodato fc NataraiarJ [200^ . here- 
after LN), we have discussed a model in which super- 
massive black holes form and grow in the center of dark 
matter haloes directly from the gas content of the halo, 
which is channelled to the center by the torques aris- 
ing from gravitational instabilities that develop in as- 
sembling massive pre-galactic discs. Models similar to 
ours, but using different sets of arguments were also pro- 
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posed in the past (lEisenstein fc Loebll9"95l: Haehnelt fc Reej 



I1993I : IKoushiappas et "ai! I I2OO4I : iBegelman et al.1 120061 ). In 



lapse of a massive quasi-star (|Begelman et al.ll2006l ). could 



par ticular, our model i mprov es upon a similar proposal 
by IKoushiappas et al.l (|2004 ) (hereafter KBD) by self- 
consistently considering the evolution of the disc surface 
density during the accretion process and by taking into ac- 
count the possibility that a very massive disc might undergo 
fragmentation and form stars rather than exclusively trans- 
port gas to the center. In this way, we obtain a robust deter- 
mination of the amount of mass accreted in the center of the 
pre-galactic disc, which is remarkably independent of many 
parameters, like the angular momentum distribution within 
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the halo, the details of the viscosity (as long as it is driven 
by gravitational instabilities) and the viscous timescale of 
the disc. As described further below, the accreted mass only 
depends on halo structural parameters (like its mass M, the 
spin parameter A and the fraction of baryonic mass that 
collapses to a disc ma), and on some well constrained gas 
dynamical parameters, like the critical value of the Toomre 
parameter for the onset of instability Qc ~ 1 — 2, and the 
gas temperature Tgas (~ 4000K for atomic hydrogen cooling, 
and ~ 300K for molecular hydrogen cooling). 

In LN we were mostly concerned with the details of 
the hydrodynamics of the disc. Here we further develop our 
model by calculating the resulting mass function of such 
central mass concentrations, available for the growth of seed 
black holes, and their volume density as a function of red- 
shift, with_thfi_standard prescription for the halo mass func- 
tion (ISheth fc Torm cn 1999). This will allow us to further 
clarify the role of the different ingredients (accretion, frag- 
mentation) of the model in the determination of the black 
hole population. The paper is organized as follows: in section 

2 we summarize the main features of our model, in section 

3 we calculate the mass function and density as a function 
of redshift, in section 4 we discuss our results and potential 
observational signatures. 



2 THE MODEL 

In the model described in greater detail in an earlier pa- 
per (LN), we consider a dark matter halo of mass M and 
virial temperature Tvir (for simplicity, the density profile 
is assumed to be that of an isothermal sphere), contain- 
ing gas which we assume to be of primordial composition, 
i.e. not enriched by metals, for which the cooling function 
is dominated by hydrogen. In the following, when numbers 
are needed, we will assume that the fraction of mass that 
collapses and fo rms a disc is Afgas ~ rridM, with rrid ~ 0.05 
||Mo et al.lll99i ). 

The angular momentum of the dark matter halo J is ex- 
pressed in terms of its spin parameter A — J\E\^^^ /GM^'''^ , 
where E is its total energy. The probability distribution of 
the spin parameter of dark matter hal oes can be obtained 
from cosmological N-body simulations (|Warren et al.|[r992l ) 
and is well described by; 
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where A = 0.05 and ax = 0.5. We further assume that the 
angular momentum of the baryonic component is a fraction 
jd ~ rnd of the halo angular momentum J. 

If the virial temperature of the halo Tvir > Tgas, the 
gas collapses and forms a rotationally supported disc, with 
circular velocity Vh, determined by the gravitational field of 
the halo. For low values of the spin parameter A the result- 
ing disc can be compact and dense and is subject to gravita- 
tional instabilities. This occurs when the stability parameter 
Q: 

CsVh ^ „ 



= 



= V2- 



(2) 



TvGT, ttGER ~ ^ ' 

where Cs is the sound speed, k — \/2Vh/R is the epicyclic 
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Figure 1. Mass available for the growth of a seed black hole as a 
function of halo mass, as obtained from Equation ||3}. The plots 
refer to the following choice of parameters: Qc = 2, Tgas = 4000K, 
md = id = 0.05, A = 0.01 (solid line), A = 0.015 (long-dashed 
line), A = 0.02 (short-dashed line). The red curves shows the 
threshold for fragmentation from Equation with ct^ = 0.06. 
Haloes on the right of the red line give rise to fragmenting discs 
in our model. 



frequency, R is the cylindrical radial coordinate, and E is the 
surface mass density and where the critical value Qc « 1 — 2. 

The subsequent evolution of the gravitationally un- 
stable disc is determined by the competition between two 
physical processes: accretion and fragmentation. Indeed, 
gravitational instabilities lead to the development of non- 
axisymmetric spiral structures (or a bar), which lead to an 
effective redistribution of angular momentum, thus possibly 
feeding a growing seed black hole in the center. This process 
stops when the amount of mass transported in the center, 
that we refer to suggestively as Mbh, is enough to make the 
disc marginally stable. This can be computed easily from 
the stability criterion in Eq. and from the disc proper- 
ties, determ ined from the halo mass and angular momentum 
()Mo et al.l il998). In this way we obtain that the mass accu- 
mulated in the center of the halo is given by: 
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where Amax = mdQc/8(md/jd)(rvir/rgas)^'^ is the maxi- 
mum halo spin parameter for which the disc is gravitation- 
ally unstable. 

It is worth noting here that the estimates obtained using 
Eq. (|3} are in fact upper limits to the actual black hole 
mass, since a fraction of the infalling gas mass might not 
be accepted by the growing black hole, if, for example, the 
rate at which the mass is delivered is significantly super- 
Eddington. This possibility is discussed further below. 

However, for large halo mass, the internal torques 
needed to redistribute the excess baryonic mass become too 



Seed black holes 3 



large to be sustained by the disc, which then undergoes frag- 
mentation. This occurs when the virial temperature exceeds 
a critical value Tmax, given by: 
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( 4qc 
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where « 0.06 is a dimensionless parameter measuring the 
critical gravitatio nal torque above which the disc fragments 
(jRice et al.ll2005l ) (note that this value has been determined 
for thin, Keplerian discs, and might be slightly different in 
the case considered here, see also LN). It is not easy to esti- 
mate how much mass can still be accreted once the fragmen- 
tation threshold is overcome. The fragmentation process is 
fast, occuring on the local dynamical timescale, and will stop 
once enough mass has been converted into stars to make the 
disc stable (i.e. the mass converted into stars will be of the 
order of Mbh, neglecting possible feedback effects), hence 
preventing further accretion, in the absence of other angular 
momentum removal mechanism. A conservative assumption 
is that, when fragmentation occurs, no mass is transported 
to the center. In the following we will discuss the resulting 
mass function for the two extreme cases where fragmenta- 
tion completely quenches accretion and where it is not taken 
into account. 

Fig. [T]illustrates the relationship between halo mass and 
potential black hole mass in our model, based on Equations 
© and dl). 

One of the main ingredients in our model is the gas 
temperature Tgas. This is computed in LN by assuming that 
the gas is in thermal equilibrium where radiative losses are 
balanced by the heating provided by the accretion process. 
We find that if the gas is atomic hydrogen Tgas ~ 4000K. 
If the gas is instead molecular hydrogen, the disc is much 
colder with « 300K. 



3 RESULTS 

3.1 The mass function of seed black holes 

From the model outlined in the previous section, using the 
mass function of dark matter haloes computed in the LCDM 
model of structure formation, we can easily obtain the mass 
function of the central mass concentration at high redshifts. 
Indeed, given the halo mass M, the probability of accumu- 
lating at the center of the disc a mass between Mbh and 
Mbh + dMBH is given by: 



P(AfBH)dMBH =p(A(M,MBH))dA, 



(5) 



where A(M, A^bh) is obtained inverting Eq. ((3)|. To obtain 
the mass function at a given redshift we need to integrate 
the probabihty given in Eq. (O over the halo mass function 
at that redshift, An/AM: 
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where the integration limits are obtained from the con- 
straints that the halo virial temperature has to be larger 
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Figure 2. Mass function of high redshift black holes. The black 
solid line refers to the following choice of parameters: z = 10, 
Tgas = 4000K and Qc = 2. The short-dashed line illustrates the 
effect of removing the threshold for disc fragmentation. The long- 
dashed line shows the effect of decreasing th critical value of Q 
to 1.5. The red line shows the results for Qc = 2 at z = 20. The 
two remaining lines refer to the cold case where Tgas = 300K at 
z = 10 (blue line) and at z = 20 (green line). 




Figure 3. Integrated mass density of black holes as a function 
of redshift. The three solid lines refer to different choices of Qc 



for the 'hot' Tg; 
and Qc ■ 
Qc = 1. 



4000K model: Qc = 3 (black), Qc = 2 (red) 
1 (green). The dashed line refers to Tgas = 300K and 



than Trr 



Tgas (in order for the gas to collapse and form 
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a disc) and has to be smaller than Tmax, obtained from Eq. 
Q, in order to prevent the disc from fragmenting. 

For the input halo mass function we adopt the stan- 
dard distribution based on the Press-Schechter formalism 
l|Sheth fc Tormen|[l999l ') and a LCDM cosmology with the 
following characteristic parameters: Qm = 0.24, = 0.76 
and as = 0.74. The resulting mass functions are plotted in 
Fig. [2] The solid line refers to the case where the gas tem- 
perature is Tgas = 4000K, the critical value of Q is Qc = 2, 
computed at redshift z — 10. The short-dashed line illus- 
trates the effects of not including the fragmentation crite- 
rion, i.e. computed for the same parameters but by removing 
the upper limit on the halo mass in the integration in Eq. 
([6|. It can be seen that fragmentation removes the highest 
masses from the population and gives an effective cut-off 
at a mass « 2 — 3 IO^Mq. The long-dashed line illustrates 
the effect of decreasing the critical value of the stability pa- 
rameter to Qc = 1.5. In this case the discs are more stable 
against gravitational instabilities and the effect is simply to 
reduce the overall normalization of the mass function. The 
red line shows the result for the same parameters as the 
black one, but at a higher redshift z — 20. As expected, at 
higher redshifts the typical masses obtained from our model 
are scaled down due to lower value of the typical host dark 
matter halo mass. The green and the blue lines refer to the 
case where the gas is assumed to cool via molecular hydrogen 
to Tgas = 300K, at z — 20 (green line) and at ^ = 10 (blue 
line). As already noted in LN, for the cold case the typical 
central concentrations have masses of the order of 10^ M© 
and are more abundant, reflecting the larger abundance of 
low mass haloes. 



3.2 Integrated BH mass density 

The integrated mass density of the central concentration as 
a function of redshift can be easily obtained from the mass 
function: 

PBh(2) = / TTT— (MBH,z)MBHdMBH. (7) 

J qMbh 

The results are shown in Fig. |3] The three solid lines refer 
to the 'hot' model with cooling dominated by atomic hy- 
drogen and Tgas = 4000K. The value of the critical Q for 
gravitational instabilities is Qc = 3 (black line), Qc = 2 
(red line), and Qc ~ 1 (green line). As expected, decreas- 
ing the value of Qc makes the discs relatively more stable 
so that the corresponding black hole density decreases. The 
dashed line refers to the 'cold' case, dominated by molec- 
ular hydrogen cooling and with Tgas = 300K and Qc = 1. 
Interestingly, the integrated black hole density in this case 
(where smaller central mass concentrations are produced) is 
comparable to the 'hot' case. This is due to the fact that in 
the cold case the population is dominated by a large number 
of small (Mbh ~ 10'' Mq) objects in the relatively more nu- 
merous small haloes, which are absent from the 'hot' case, 
because the gas in such small haloes has a temperature larger 
than the virial temperature and does not collapse to form a 
disc. 

The present day black hole density, as derived 
by Soltan's argumen t is between 2 — 4 lO^Mo/Mpc^ 
jYu fc Tremainj|2002l ). It would then naively seem that 
our model might overestimate this quantity at low red- 
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Figure 4. The relationship between spin parameter and black 
hole mass for a halo of mass M = 3 IO^Mq at 2 = 14, according 
to the LN model (solid line) and to the KBD model (dashed line) . 
For the KBD model, the additional parameter tvisc was assumed 
to be lOMyrs, while for the LN model we have assumed Qc = 2. 

shifts. It should be however kept in mind that our proposed 
model strictly offers a quantitative estimate of the black hole 
mass function only at high redshift, before the intergalactic 
medium is enriched by metals. It is then interesting to see 
that our model is able to provide a reasonable black hole den- 
sity of a few times Mq /Mpc^ at a redshift, say, 2 = 15 to 
allow for further AGN-like growth to the present day density. 
The consequences of using the seed mass function derived 
above to compute the supermassive black hole mass function 
at low redshift using the merger-tree approach is currently 
underway (Volonteri, Lodato & Natarajan, in preparation). 

Another important aspect to be taken into considera- 
tion is when the transition occurs from mostly molecular gas 
to mostly atomic gas. It is reasonable to assume that at very 
high redshift the gas is predominantly molecular becoming 
atomic after the birth of the first st ars, when the UV ba ck- 
ground has increased substantially (|Oh fc Haiman|[2003l ) . A 
detailed investigation of this problem is however, beyond the 
scope of the present paper. 



4 COMPARISON TO EARLIER WORK 

A model similar to our own has been also discussed by KBD. 
A detailed discussion of similarities and differences between 
the two models is provided in LN. Here we point out how 
these differences affect the resulting mass function. There 
are two main differences between the two models: (i) we 
include the evolution of the disc surface density due to the 
redistribution of angular momentum caused by gravitational 
instabilities, while KBD only consider the evolution induced 
by the infall of matter onto the disc, (ii) We include the ef- 
fects of disc fragmentation, which are neglected by KBD. As 
a result of point (i) above, the KBD results depend on some 
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additional free parameters which are implicitly (and self- 
consistently) included in our model, like the viscous time- 
scale in the disc. These affect the functional relationship 
between central mass and spin parameter for a given halo 
meiss dramatically. In Fig.|4]we plot Mbh as a function of A 
for a given halo mass according to our model and the cor- 
responding plot obtained using the KBD model (see Figure 
caption for details). We note that for the same structural 
parameters our model predicts lower masses than KBD. 

The relationship between Mbh and A according to KBD 
is a simple power-law, so that for a given halo mass, the 
black hole mass distribution is essentially log-normal, as the 
spin parameter distribution, while this is not the case in our 
model. Our results suggest that dark matter haloes with 
high spin values (above a certain value of the spin parame- 
ter A) do not host black holes at early times. This is because 
for a given M, if the disc is not compact enough (or con- 
versely, for a given A, if the disc is not massive enough) it is 
stable to gravitational instabilities. Such haloes likely have 
more conventional seeds from remnants of the first genera- 
tion of stars or eventually host black holes as a consequence 
of mergers at later times (Narayanan ct al. 200^). 

Additionally, the black hole mass function as derived 
from the LN model is also skewed to high black hole masses, 
due to the effective mass cut-off set by the onset of fragmen- 
tation in very massive discs, as shown in Fig. [2] 



5 CONCLUSIONS AND DISCUSSION 

One of the key assumptions made in this work is that the 
centrally concentrated mass in pre-galactic discs eventually 
forms a black hole. As mentioned above in Section 2, not 
all of the mass that flows into the center of the halo might 
be accepted by a growing black hole, so that our mass esti- 
mates should, strictly speaking, be considered upper limits. 
We wish to stress that in this paper we are mostly con- 
cerned with determining the total mass budget of such sys- 
tems, rather than in the eventual fate of the accreting mat- 
ter. Detailed models of the ultimate fate of the gas as it 
forms a seed black hole have shown how t h e growth can be 
super - Eddington (|Volonteri fc Reej l2005l : iBegelman et al.l 
I2OO6I ). Alternatively, the accretion flow at the smallest scale 
might bre ak down and become chaotic. In such a situation 
l|King fc Pringle 2(M ) the accretion process does not spin 
up the hole efficiently and th e accretion efficiency is thus 
kept low, to e ~ 0.05 (see also lLodato fc Pringle|[2006l ). ef- 
fectively increasing the Eddington accretion rate. 

Once a seed black hole with a mass of a few times 10* — 
10^ Mq has formed (cf. the mass function shown in Fig. [2] at 
z — 20), it can initially keep accreting following our model. 
A strict upper limit to the accretion rates through the disc 
in our model is given by the threshold for fragmentation, 
which is ~ O.OlMo/yr (see LN), which corresponds to the 
Eddington li mit for a 10^ black hole, with an efficiency 
of e = 0.05 (|King fc Pringld bOOd ). Therefore, we can be 
confident that in the majority of cases most of the available 
mass will indeed be accepted by the growing black hole, so 
that our upper limits for the black hole density are fairly 
realistic estimates of the actual values. 

Our simplified analytical model for accretion essentially 
considers isolated haloes, while the Press-Shechter formal- 



ism adopted to determine the mass function is based on the 
merger history of the dark matter haloes. Mergers can have 
two effects on the results described here. Firstly, a merger 
might disrupt the disc before it has enough time to form and 
accrete to the center. These two processes occur on the same 
timescale, i.e. the internal dynamical time of the halo, which 
is generally shorter than the survival time of the halo (see 
iLacev fc Colli 1991 Fig. 5), so our model is consistent from 
this point of view. A second point is that a halo at a given 
redshift might already contain a hole formed previously in 
the parent haloes. However, as long as the growth of the seed 
is dominated by accretion rather than by merger of different 
seeds (which is generally the case, see IVolonteri et~al]|2003l ) 
we can still use our formalism to determine the new seed 
mass. At a given redshift, some seed black holes will have 
formed anew from the lowest mass haloes, while some oth- 
ers (contained in the most massive haloes) will already have 
formed earlier but they will still have grown to a new, much 
larger mass, that can be computed from Eq. ([3J. Our model 
cannot distinguish between the two cases, and our mass 
function includes both contributions. In order to disentangle 
the two one would need to follow the detailed merger history 
of individual haloes rather than simply using the halo mass 
function at a given redshift as an input (an approach we 
pursue in Volonteri, Lodato fc Natarajan, in preparation). 

Massive black holes may also be related to the re- 
cent discovery of ultra-lu minous x-ray sou r ces (ULXs) 
in n e arby galactic discs (|Mushotzkvl |2004| : iMiller et al.l 
12004 iDewangan et al.ll2006l ). The observed ULX luminosi- 
ties of ~ lO'^^ergs"^, are inconsistent with s tellar mass 
black holes unless they a re highly beamed (Ki ng et al.l 
I2OOII : iBegel man et al]l2006l ). Intermediate mass black holes 
(IMBHs) in the mass range 10^ M© < A/bh < 10* M©, 
have the appropriate Eddington luminosities to explain 
ULXs. While stellar collapse is unlikely to pro duce this mass 
range of black holes (|Frver fc Kalogerall200ll) . our proposed 
mechanism provides a natural way to do so. Recent detailed 
numerical studies suggest that these IMBHs mi ght be criti- 
cal to the build-up of sup ermassive black holes (IMicic et al.l 
2006; Fro geau et al.ll2006l '). A fraction of IMBHs are expected 
to be wandering in the haloes of galaxies as a result of fre- 
quent merging. If any of these encounter gas-rich regions 
they are likely to accrete and be detected as ULX's. 

Accretion onto the seeds produced via the disc accretion 
scenario described a bove at extremely hig h redshifts could 
power mmi-quasars l|Madau fc ReesH200ll ) that likely play 
a significant role in the reionization of the Universe. These 
sources provide an attractive way to partially reionize the 
low density inter-galactic medium, as the harder radiation 
emitted from these mini-quasars in the X-ray is more likely 
to escape compared to the softer UV radiation from mas- 
sive stars. Miniquasars are ~ 10 — 100 times more efficient 
at producing ionizing radiation per processed baryon than 
metal-free Pop-Ill stars, therefore require a much smaller 
reservoir of cold gas to reionize the Universe. 
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